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METHOD OF PROCESSING OSCILLATORY DATA 

The present invention relates to a method of processing oscillatory data. It is 
particularly, but not exclusively, concerned with a method of processing oscillatory 
response data, and in particular non-stationary responses caused for example by a 
characteristic frequency of a system changing with time. 

There are many situations where the response of a system, such as a gas 
turbine engine or a component thereof, which is subjected to random or engine order 
(EG) excitations, displays non-stationary characteristics. These may occur even where 
the excitation itself is reasonably stationary. 

For EG excitation, engine speed changes will cause a varying response level as 
the EG frequency approaches or moves away from a modal frequency. Thus the 
response of the system will be transient in nature and non-stationary in that sense. If, 
in addition, the modal frequency in question is itself changing with time or engine 
speed, then the transient response will be further affected, resulting in a complex 
speed, or time, response profile. 

Random excitation data is normally considered stationary if statistical properties 
such as long time averaged root-mean-square level in a frequency band are constant 
with time. The non-stationary case can arise if the modal frequencies are changing 
with time. 

Changes in modal frequency can occur for a variety of reasons. For exanriple, in 
a gas turbine engine, these variations can be due to thermal and/or centripetal force 
stiffening. Within missile systems these variations may also be due to mass loss 
effects. 

In mechanical systems, damping provides a measure of the energy dissipation 
capacity of a structure or system. The measurement of damping can be useful 
because it aids the understanding of the behaviour of a structure, and knowledge about 
the level of damping can be valuable in troubleshooting and assessment of potential 
and actual problem fixes. Damping can also be used to provide parameter values for 
and validation of models. 
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Changes in damping over time may be useful to indicate the "health" of a 
system, and can also be used to provide an indication of potential instability of a system 
such as "flutter". 

It is difficult to determine the characteristics of a mode exhibiting such non- 
stationary behaviour, arising due to modal frequencies changing with time. For 
example it is difficult to determine the modal bandwidth or damping associated with the 
mode. 

This is because conventional analysis methods and systems do not adequately 
handle non-stationary data such as that discussed above, and therefore it is difficult to 
obtain accurate measurements of any of the above characteristics in systems 
displaying non-stationary characteristics. 

In practice the most common approach used for dealing with non-stationary 
random excitation data at the moment is the segmented approach outlined in Random 
Data: Analysis and Measurement Process by J.S. Bendat andA.G. Piersol; Wiley 
Interscience, 1971. We are not aware of any conventional approach which deals with 
EO data and compensates for changes in the non-stationary mode characteristic data. 

The segmented approach splits the oscillatory response data into segments, 
each of which are considered to be stationary and therefore processed as such in the 
normal manner. The above reference describes why results from this approach should 
be treated with caution and may only be useful in a qualitative sense. 

Specifically, a problem arises with the short: time interval bias error arising from 
changes in the characteristics which occur within the period of an individual segment, 
e.g. if the modal frequency changes significantly relative to the modal bandwidth within 
a segment. In order to suppress the time interval bias error, a short averaging time 
(subrecord length) T is required, but to obtain the desired spectral decomposition, a 
narrow resolution bandwidth B© is required. The result is a relatively small Be.T product 
and hence a large random error. 

However, in many cases the requirement for adequate frequency resolution will 
mean that segment durations are long, with significant modal frequency changes 
relative to the modal bandwidth and therefore a serious bias occurs in the bandwidth or 
damping estimates. 
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Therefore the segmented approach is in principle problennatic and in some 
situations may not even allow qualitative results to be obtained. 

The present invention seeks to address some or all of the problems associated 
with the segmented approach. An aspect of the present invention also seeks to provide 
a methodology for handling EO data. 

In general terms, the present invention provides a method which permits the 
transformation of non-stationary response measurements into stationary data for 
analysis. 

The present invention relates to oscillatory systems in general, but will be 
described and illustrated for oscillatory data which specifically describes vibration 
phenomena. Alternative applications of the invention, for example in electrical networks 
in which it is desired to determine the bandwidth when the characteristic frequency of 
the circuit is changing with time, are equally possible. 

A first aspect of the present invention provides a method of processing 
oscillatory response data from a resonant system comprising: 

obtaining data measuring an oscillatory response of the system; 

estimating the variation in natural frequency of a mode of said response; 

filtering the data around a selected frequency to obtain a filtered response; 

determining a carrier signal whose frequency variation with respect to time is 
equal in magnitude to said estimated variation in natural frequency; and 

modulating the amplitude of said carrier signal using said filtered response to 
obtain a modulated carrier signal. 

The selected frequency around which the data is filtered may be the natural 
frequency of a mode of the response (and preferably the mode which is the subject of 
the estimating step) or the frequency of an engine order. In frequency terms, the 
filtering may be static (e.g. centred on a fixed frequency) or dynamic (e.g. centred on a 
varying frequency). 

The frequency variation of the carrier signal may be equal to the estimated 
variation in natural frequency or may be the negative of the estimated variation. 

Preferably the frequency of the carrier signal is greater than the difference 
between the highest and lowest values of the natural frequency of said mode over the 
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period of interest. This may allow the data from two sidebands formed in said 
modulated carrier signal to be analysed separately with no overlap between them. 

Transforming the typically non-stationary oscillation data according to the 
method of the present aspect may make it effectively stationary, and thus in a more 
useful form for further analysis. 

The step of estimating the change in natural frequency may include obtaining 
multiple Fourier transforms of the measured data, for example as described in "The 
automatic interpretation of vibration data from gas turbines", Allwood, King & Pitts, The 
Aeronautical Journal of the Royal Aeronautical Society, March 1996. 

Alternatively to the methods described in this reference, the same zmod 
(waterfall or Campbell diagram) data containing the multiple Fourier transform data may 
be analysed to determine likely natural frequency values via a standard modal curve 
fitting approach and/or by manually cursoring likely estimates prior to conducting a least 
squares type fit to the data to determine a smooth time frequency profile. 

Alternatively, for example where the mode is remote from others and the 
response may be bandpass filtered around the mode, the step of estimating the change 
in natural frequency may include calculating a running average of the instantaneous 
frequency of the response. 

Furthermore, model data or a combination of model data with experimental data 
may be used in the step of estimating the change in natural frequency. 

The step of estimating the variation in the natural frequency of the mode of the 
response may also include further steps which fine-tune the estimation. For example, 
an optimisation approach may be used to alter the natural frequency variation, under 
predetermined constraints, in order to minimise the modal bandwidth associated with 
the transformed data and mode of interest. The constraints may impose limits on the 
smoothness associated with the functional form of the natural frequency estimate 
versus time. For example, changes which are considered not to be smooth enough, or 
which are considered to be too quick might be rejected as unrealistic. 

When the data being processed comes from engines or, in particular, gas 
turbines, in many cases it will have arisen from conducting engine speed-up or slow- 
down tests in which the engine speed changes over a given time interval. Under these 
conditions it is common for the natural frequency variation with time to be approximated 
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by a quadratic equation because of the long time constant associated with the thermal 
or centrifugal-force-stiffening effects causing the natural frequency changes. Higher 
degree polynomial approximations or other functional forms or numerically described 
profiles could also be used. Numerically described profiles could, for example, be 
5 obtained by tracking the peak on a running spectrum. 

An advantage of the method of processing of this aspect is that modes which are 
close in frequency terms to the mode being identified are effectively "smeared" In 
frequency by the time averaging process if their natural frequency time profile differs 
from that of the mode being identified. This increases the likelihood of the mode of 

10 interest being identified, and allows better estimation of modal bandwidths from a single 
mode estimation methodology. 

Another aspect of the present invention provides a method of analysing a 
resonant system, including the steps of: processing according to the above aspect; and 
analysing the modulated carrier signal to determine one or more characteristics of the 

15 system. 

As a result of the transformation of the response data, standard analysis 
techniques may be used on the modulated carrier signal to determine features of 
interest such as bandwidth and damping, and when the data being analysed is 
effectively stationary, the analysis may be done with high frequency resolution and long 
20 averaging times, without bias errors resulting from natural frequency changes. 

For example, the step of analysing may include determining a power spectral 
density (PSD) function of the modulated carrier signal. 

Alternatively, if the input, which is driving the system, is known, a frequency 
response function computation may be conducted. Again due to the processing of the 
2 5 data, standard random data analysis techniques may be used. 

A further advantage of the methods of the present invention is that bias (or other 
systematic) errors in the natural frequency estimation may be tolerated since the modal 
bandwidth is unaffected by the transformation and bias errors in the natural frequency 
simply affect the apparent frequency of a mode which has been made stationary. 
30 For EO excitation a further advantage of the methods of the present invention is 

that it enables a modal curve fitting approach to be used on the transformed data, even 
if the EO does not fully traverse the mode; i.e. the skirts or only part of the resonance 
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curve profile may be analysed to determine damping, even when the engine speed falls 
short of causing the relevant EO frequency to equal and exceed the modal frequency. 

The methods of the previous aspects may conveniently be implemented in 
software, for execution on any appropriate digital computer. The software may also 
embody preferred or optional features of the methods of the previous aspects. The 
methods may be performed on-line, or off-line on stored measurement data. 

Thus further asfjects of the invention respectively provide a computer system 
operatively configured to implement any of the methods of the previous aspects of the 
invention; computer programming product or products (such as ROM, RAM, floppy 
discs, hard drives, optical compact discs, magnetic tapes, and other computer-readable 
media) carrying computer code for implementing any of the methods of the previous 
aspects of the invention; and a computer program perse for implementing any of the 
methods of the previous aspects of the invention. 

Further aspects of the present invention also provide apparatuses for carrying 
out the methods of the above aspects. 

In particular, these aspects provide an apparatus for processing oscillatory 
response data from a resonant system, the apparatus including: 

a sensor for measuring an oscillatory response of the system; and 
a processor which is adapted to: 

receive measurement data from the sensor; 

estimate from the data the variation in natural frequency of a mode of said 
response; 

filter the data around a selected frequency to obtain a filtered response; 

determine a carrier signal whose frequency variation with respect to time is equal 
in magnitude to said estimated change in natural frequency; and 

modulate the amplitude of said carrier signal using said filtered data. 

In an apparatus for analysing the resonant system, the processor may be further 
adapted to analyse the modulated carrier signal to determine one or more 
characteristics of the system. 

The processor may be implemented in dedicated hardware, or it may be a 
suitably programmed computer. 
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Preferred or optional features of the methods of the previous aspects may be 
embodied in corresponding preferred or optional features of the apparatus of this 
aspect. 

The system of any of the above aspects may be a model system. Alternatively 
the system may be a mechanical system, such as a gas turbine engine or a component 
thereof. Alternatively the system may be an electrical network. 

An embodiment of the present invention will now be described in relation to the 
accompanying drawings, in which: 

Figure 1 shows an oscillatory response of a model system; 

Figure 2 shows a power spectral density (PSD) computation of the response of 
Figure 1; 

Figure 3 shows the response of Figure 1 which has been band-pass filtered 
around a mode of interest; 

Figure 4 shows an expanded part of the response of Figure 3; 

Figures 5a-5c show schematic representations of, respectively, a magnification 
of the narrowband response, a data set of signal zero crossing separation times of that 
response, and a data set of instantaneous frequencies of that response; 

Figure 6 shows the instantaneous frequency of the response over the period of 
analysis; 

Figure 7 shows the variation in the carrier frequency used in the method of a first 
embodiment of the present invention; 

Figure 8 shows a PSD calculation of the response after the method of an 
embodiment of the present invention has been used; 

Figure 9 shows an oscillatory response of a second model system; 

Figure 1 0 shows a zmod or waterfall plot of the response of Figure 9; 

Figure 1 1 shows an estimate of the natural frequency of the mode to be 
identified from the response of Figures 9 and 10 over the period of analysis; 

Figure 12 shows the variation in the carrier frequency used in the method of a 
second embodiment of the present invention; 

Figure 13 shows the variation in the carrier phase used in the method of the 
second embodiment; 
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Figures 14 and 15 show the carrier signal used in the second embodiment over 
an interval of 0.04 seconds at respectively the start and end of the simulation period; 

Figures 16-19 are expanded plots of the time interval of 50.0 to 50.3 seconds in 
the second embodiment, showing respectively: the response; the bandpass filtered 
response; the carrier signal; and the resulting modulated signal; 

Figure 20 shows the zmod plot of the resulting modulated signal produced by the 
second embodiment; 

Figure 21 shows a PSD calculation of the response after the method of the 
second embodiment has been used; 

Figures 22 and 23 respectively show the magnitude and phase of the frequency 
response function resulting from the processing of the second embodiment; 

Figure 24 shows schematically the zmod or waterfall plot of the response data 
used in a third embodiment of the present invention; 

Figure 25 shows the forced response amplitude profile for the simulation used 
for the third embodiment; 

Figure 26 shows the response amplitude profile after a method according to the 
third embodiment of the present invention has been used; 

Figure 27 shows the frequency variation of the response of the upper sideband 
of the response after processing according to the third embodiment; and 

Figures 28 and 29 show schematic zmod plots of data respectively prior to and 
after processing in order to illustrate an aspect of the theory. 

General Theory 

A general theory of embodiments of the present invention is set out below. 

For this purpose, the systems to which the embodiments of the invention are 
applied are considered to have an input x(t) and an output y(t), both of which are 
functions of time t. The output may represent a measured or predicted response. 

A mode of interest of the system has an angular modal frequency, the natural 
frequency, which also changes with time. An approximate description of its angular 
frequency profile (jOn(t) is obtained, for example by processing y(t) by zmod plotting, 
bandpass filtering, zero crossing analysis, etc., and/or making use of model data. 

It is desired to transform y(t) to compensate for the non-stationary behaviour 
associated with the frequency changes of the mode. The invention accomplishes this 
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by amplitude modulating y(t) witii a carrier signal whose frequency also changes, in a 
synchronous manner with the mode frequency. Although the term carrier signal is used 
in this specification, the signal formed is quite different from normal carrier signals, 
formed for example in communications technology, in that it varies with the natural 
frequency. 

Thus a varying carrier angular frequency is created: 
coc(t) = cOs - [a)n(t) - a)n(t*)] Equation (1) 

where: C0n(t*) is the value at some arbitrarily chosen instant (e.g. t* = 0 or when 
the response peaks), and 

cos is a constant shift frequency that places signal components associated with 
the mode of interest of the transformed signal into another part of the spectrum, cus is 
chosen so as to be sufficiently large to ensure that the lower and upper sideband 
components of the modulated response (see below) are adequately separated in 
frequency. 

Note that since cos and a)n(t*) are both constants, they could in fact be 
incorporated into a single constant term (a)constant); ©-g. the carrier frequency then 
becomes 

a)c(t) = Wconstant " Cjan(t). In somo applications it may be more convenient to think in terms 
of the form expressed in equation (1) although both are equally valid. 

Similarly, the carrier frequency for the lower sideband formulation (described 
below) could equally be defined as cjOc(t) = ooconstant + Wn(t). 

The phase 0c(t) of the carrier signal is obtained by integrating a)c(t). The phase 
at t = 0 may be arbitrarily set (e.g. to zero). The integration may be done by a variety of 
well known techniques, such as the first order Taylor approximation for a discrete 
process: 

ec(t) « Gc(t - dt) + {dt. C0c(t - dt)} 

where dt is the time interval of integration. 

The carrier signal then becomes: 

Sc(t) = sin {ec(t)} 

although of course a cosine function could also have been chosen since the 
initial phase is arbitrary - the requirement being simply that the rate of change of the 
phase Gc(t) equals the carrier frequency (jOc(t) as defined above. 
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A modulated or transformed signal ym(t) is then obtained by multiplying the signal 

y(t) by Sc(t): 

y^(t) = sc(t).y(t). 

This results in a time dependent transformation of frequencies in a continuous 
manner, enabling further analysis using standard stationary data analysis techniques to 
yield bandwidth and damping information. The mode of interest now appears to be 
stationary with an angular frequency equal to [oos + Wn(t*)] when looking at the upper 
sideband data of the modulated signal. 

The basic modulation effect and the "upper" and "lower sideband" terminology 
can be appreciated by considering the simple case where Sc(t) = sin (0) and y(t) = cos 
(cp), in which case ym(t) = Sc(t).y(t) = y2.[sin (9-(p) + sin (e+cp)]. (G-cp) is referred to as the 
lower sideband component and (9+(p) is referred to as the upper sideband component 
of the modulated signal. The upper sideband component has an angular frequency 
which is the time rate of change of (9+(p), which in turn is equal to the time rate of 
change of each of 0 and cp separately - i.e. the sum of the instantaneous frequencies 
associated with the carrier and original response signals. 

Thus, for the case where the signal y(t) has a time varying frequency component 
of {u}n(t) + A} where A is constant, then (from equation 1) and the above it can be seen 
that this component is transformed to an upper sideband frequency of: 

{a)n(t) + A} + {uJs - [con(t) - Wn(t*)]} = {Ws + C0n(t*) + A}, which is a constant, as it 
has no time varying components. 

Therefore, frequency components equal to the natural frequency {A = 0} are 
transformed to the constant frequency {Ws + CA)n(t*)}, and components with a given 
(constant) separation frequency from the natural frequency are transformed in a 
manner that maintains that separation frequency. 

The formulation of Equation (1) above is thus referred to as the "upper sideband 
formulation" in that the upper sideband created in the modulated signal is stationary for 
the purposes of identifying the mode of interest. Using a carrier frequency described 
by: 

Wc(t) = Ws + [a)n(t) - Wn(t*)] 

would generate a modulated signal whose lower sideband would have stationary 
properties for the mode of interest (the "lower sideband formulation"). The first 
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embodiment described below uses the lower sideband formulation whilst the second 
and third embodiments use the upper sideband formulation. 

in situations where y(t) is not sufficiently narrowband, it can be bandpass filtered 
prior to modulation to prevent interference between signal components which may be 
transformed into the same frequency range; e.g. a low frequency component fl in y(t) 
may be translated up to f2 in the modulated signal but this may coincide with a high 
frequency component f3 in y(t) which has been translated down and now also appears 
as f2 in the modulated signal. 

Therefore in most cases, the modulated signal will be determined as: 

ym(t) = Sc(t).y(FL, FU, t), 

where y(FL, FU, t) represents the bandpass filtered result after filtering y(t), and 
FL and FU respectively denote the lower and upper passband frequencies and may be 
constant or time dependent. FL and FU are generally selected to ensure that frequency 
components near the mode of interest are unaffected whilst frequencies away from this 
local region are rejected. 

The selected sideband of the modulated signal can now be further processed 
using standard techniques for stationary data, such as Power Spectral Density (PSD) 
function computations, to yield information about the modal bandwidth and damping. If 
further processing is to be done in the frequency domain, then the signal components in 
the other sideband are ignored. If further processing is to be done in the time domain, 
then the modulated signal may be bandpass filtered around the components in the 
sideband of interest prior to any such processing. 

If absolute, rather than relative, magnitude values are important, the modulated 
signal can be multiplied by a factor of two to correctly scale the levels. 

In situations where the input x(t) is known, then this may also be modulated in a 
similar fashion to the response y(t). Thus a modulated input signal Xm(t) = Sc(t).x(t) is 
formed. As with the response signal y(t), the input signal x(t) may be bandpass filtered 
first to make it sufficiently narrowband. The two modulated signals (input and output) 
may then be used, for example to compute the frequency response function between 
the signals, using standard spectral and cross-spectral techniques as a precursor to 
mode and damping estimation. 
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The above processes may be carried out in an analogue or a digital manner. 
Multiple inputs and/or outputs may be processed using the above methodology. 

If the excitation is also non-stationary so that the time average force level varies 
significantly, then this may also be compensated for if required. Generally the 
excitation changes have a time constant that is much longer than any considered 
vibration time constants (i.e. the changes are slow in comparison to the vibration itself). 
One approach is to track the mode of interest and compute a smooth running root- 
mean-square (rms) profile, a(t), of the response. This may then be used to normalise 
the modulated signal prior to any standard analysis. The signal is thus made stationary 
as either the upper or lower sideband part of {[Sc(t).y(FL, FU, t)]/[a(t)]}, depending on 

the formulation used. 

In small time regions where the modal frequency of the mode of interest 
coincides with other modal frequencies, the transformed data may be ignored if 
necessary prior to performing any standard stationary data analysis. 

As the natural frequency of the non-stationary sideband will vary to a greater 
extent after processing as described above, it is possible that data from the two 
sidebands overlaps at some point in time. This is generally undesirable since such an 
overlap can affect the analysis of the stationary sideband. 

Therefore it is preferable that the carrier signal is chosen so that such an overlap 
does not occur in the period of interest. This can be achieved by choosing a carrier 
frequency which is greater than the difference between the highest and lowest values of 
the natural frequency of said mode over the period of interest, as shown below. The 
analysis below is based on positive values of the carrier signal frequency Ws, although 
negative values are possible, and the same considerations would apply, but with the 
"upper" and "lower" sidebands being reversed. 

The zmod plots of the response of a fictional system with a single mode whose 
natural frequency changes, and which is subject to random excitation are shown in 
Figures 28 and 29 which respectively show the data before and after processing 
according to an embodiment of the present invention. 

The illustrated transformation achieved by the processing of this embodiment 
can be considered a "good" transformation in the context of achieving separation of the 
upper and lower sidebands by appropriate selection of Ws, such that the upper 
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sideband (USB) frequencies are separated from the lower sideband (LSB) frequencies 
over at least a significant portion of the time history, and preferably over the entire 
period of interest. 

To determine the values of cOs which may allow such "good" transformation, we 
consider the two frequency profiles: 
profile 1 : u)i(t) = cjOn(t) + 5 
profile 2: a)2(t) = {jOn(t) - 5, 

where a)n(t) describes how the natural frequency varies, and 6 represents the 
maximum frequency separation from the natural frequency which may be of interest. 
Typically, 5 will be of the order of the modal bandwidth. 

Under the transformation of embodiments of the invention cjOi(t) gets transformed 
to cjo/(USB,t) and (jD/(LSB,t); and waW is similarly transformed. Using the superscript 
T' to represent the transformed data, we will therefore have the following profiles 
associated with the transformed data: 

profile 3: (jo/(USB,t) = coi(t) + a)c(t) 

profile 4: oo/(LSB,t) = coi(t) - oac(t) 

= con(t) + 5 - COc(t) 

(after substitution from above) 
profile 5: (ja2"^(USB,t) = C02(t) + {jDc(t) 

= CjOn(t) - 6 + CjOc(t) 

(after substitution from above) 
profile 6: a)2^(LSB,t) = a)2(t) - (jOc(t). 

where USB & LSB refer respectively to the upper and lower sidebands. 

If required, the data may be initially filtered (prior to the transformation) to ensure 

> 

that the signal only has significant frequency components in the approximate range 

{A)n(t) +/- 8. 

Imposing a requirement that the USB & LSB frequencies are separated (i.e. do 
not mix) over the entire period of interest or time history implies that: 
co2"^(USB.t) > a)/(LSB,t) for all t. 

By making all the relevant substitutions for 032^(USB,t) and (^/(LSB.t) from 
above, we obtain the-^condition that: 
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(ji)n(t) - 8 + Wc(t) > u)n(t) + 8 - Wc(t), and therefore that: 
U)c(t) > 8. 

Substituting from Equation (1) gives: 
Ws > 8 + [a)n(t) " a)n(t*)] 

Denoting the maximum and minimum values of (jOn(t) respectively by con^max and 
cjOn.min, the above preferred constraint equation on cos may be conservatively and further 
simplified to: 

cOs > 8 + [con.max " Wn(t*)] and thence to 

Ca)s > 8 + [COn.max " CA^n.min] 

which are sufficient conditions. 

In cases where the signal bandwidth (due to significant natural frequency 
changes) is much larger than the modal bandwidth, i.e. where 

[Wn,max " Wn.min] 8, 

the constraint equation approximately simplifies further to: 



Therefore selecting oos to be greater than the difference between the highest and 
lowest natural frequency for the mode of interest is sufficient to ensure that the lower 
and upper sidebands are sufficiently separated. Clearly setting (jJs > (J^n.max also 
satisfies the condition. 

Even if this most stringent condition is not met, the transformation effected by 
embodiments of the invention may still be useful as there will normally be parts of the 
data time history which reveal useful information. 

A similar construction can be done for the Engine Order situations. Again we 
consider an upper sideband formulation and use similar reasoning and notation to that 
above. 

Setting a)n(t*) = 0, then u)c(t) = cOs - cji)n(t). 

Consider the transformation of the engine order frequency (ooeo): 




a)eo"^(USB.t) 



Weo(t) + a)c(t) 
U)eo(t) + U)s - (A)n(t) 



Weo"^(LSB,t) 



Weo(t) - aJc(t) 

= a)eo(t) - COs + COn(t) 
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The condition for a)s is 
a)eo"^(USB,t) > a3eo^(LSB.t) 
Therefore we obtain: 

oJeo(t) + 00s - cjan(t) > a)eo(t) - cjOs + Wn(t) and 

Again Ws > Wn.max would of course be conservative and satisfy the condition. 

First Embodiment 

A modelled scenario is used to illustrate a method according to a first 
embodiment of the present invention, 

A model system is subjected to random vibration excitation. The model system 
has two dominant modes, both of which are time-varying, that create a typical 
response, as shown in Figure 1. In a real situation, this response could be, for 
example, a strain or acceleration measurement from a mechanical system such as a 
gas turbine engine or a component thereof. 

The simulation is conducted over a period of 100 seconds with a time interval 
between discrete points of the simulation of 1/5000 seconds. 

The two modes have the following characteristics: 

The natural frequency of Mode 1 varies from 100 Hz to 1 10 Hz between 0 and 
80 seconds and then remains constant. Mode 1 has Q = 100, Q being the quality factor 
which can be calculated as the characteristic or natural frequency divided by the modal 
bandwidth. 

The natural frequency of Mode 2 varies from 200 Hz to 250 Hz between 0 and 
80 seconds and then remains constant. Mode 2 has Q = 50. 

The -3dB modal bandwidths associated with these modes are approximately 1 
Hz and 4.5 Hz respectively. 

The system response is very clearly non-stationary as a result of the large 
changes in natural frequency. The PSD computation of the whole period of analysis in 
Figure 2 shows the smearing of the frequencies due to this non-stationary behaviour. 
The frequency resolution of this analysis is approximately 0.1562 Hz. 

To estimate the variation of the natural frequency of Mode 1 , the system 
response is preferably bandpass filtered (e.g. between 50-150 Hz) to obtain principally 
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the contribution of Mode 1 to the response. This filtered response is shown in Figure 3 
and a part of it is shown in magnified form in Figure 4. 

This filtered response may be analysed to compute the "instantaneous" 
frequency and its variation over the time of the response. This analysis may be done 
by computing the time between successive zero-point crossings, which are taken to 
represent half periods of the response. 

Figure 5a shows a schematic magnification of a narrowband response, showing 
the zero crossing points. From the zero crossing points, a time-varying data set of 
separations can be obtained, as shown in Figure 5b. This can be translated to a time- 
varying data set of instantaneous frequencies (fj = 1/2ti) as shown in Figure 5c. Note 
that the perturbations within the small time sample shown are generally not significant 
enough to be noticeable on the overall plot of the time separation or instantaneous 
frequency. 

Alternatively, local curve fitting analysis can be used to compute the 
"instantaneous" frequency and its variation. For example a sinusoidal function (in the 
interval [0, it] which describes a half cycle) may be fitted to all data points between 
consecutive zero-crossing points. This fitting is repeated for all points to give the 

frequency profile. 

The frequency data may itself be filtered. A running average of the 
instantaneous frequency or a curve fitting process may then be used to estimate the 
variation of the natural frequency of the mode over the period of the response. Figure 6 
shows a plot of the variation of the instantaneous frequency Gagged line) and of the 
smoothed frequency (obtained using curve fitting) with respect to time. 

Alternatively an estimate of the natural frequency variation could be obtained by 
analysis of a waterfall type spectral analysis of the whole response signal as shown in 
the second embodiment below. 

Next, a shift frequency of, in this case 50 Hz, is chosen. That is, a carrier signal 
is generated which has a frequency of 50 Hz more than the change in the smoothed 
instantaneous frequency of Mode 1 over the course of the response. The magnitude of 
the shift frequency ensures that the lower and upper sidebands created by the 
modulation process are well separated. 



wo 2005/068949 PCT/GB2004/005055 

- 17- 



The carrier frequency of the first embodiment is formed using the lower sideband 
formulation described above and is shown in Figure 7. 

The carrier signal thus formed is then multiplied by the narrowband filtered 
response (shown in Figure 3). 

The resulting signal can then be processed using standard analysis techniques 
and methods for stationary data sets. 

For example, the PSD of the resulting modulated signal may be computed, as 
shown in Figure 8. The lower frequency sideband indicated represents the spectral 
distribution associated with Mode 1 after the transformation to the stationary data set 
and mapping to a lower frequency. The other spectral components are the upper 
frequency sideband of Mode 1, which, have been caused to become even more non- 
stationary by the processing. 

As can be seen from a comparison of the PSD shown in Figure 8 with that 
shown in Figure 2, the signal is much less spread in frequency and therefore can be 
analysed more easily using standard methods (which need not include determining a 
PSD function). 

Specifically, the response due to mode 1 has been effectively made stationary 
without energy being distributed over a wide band due to natural frequency changes. 

In particular, features such as the vibration level and modal bandwidth and/or 
damping information of the lower frequency sideband can be much more accurately 
and easily determined. 

Second Embodiment 

A second modelled scenario is used to illustrate a method according to a second 
embodiment of the present invention. 

A two mode system is excited by random broadband noise. The time history 
response is shown in Figure 9. The running spectral profile of this response data, 
commonly known as the "zmod", waterfall or Campbell diagram, is shown in Figure 10. 
Such plots usually plot frequency on the vertical axis, but in the particular plotting 
program used to produce Figure 10, the axis represents the non-dimensional parameter 
of the frequency (f) divided by the frequency resolution (Af). 

The "zmod" or waterfall is obtained from a well known analysis technique which 
involves splitting the oscillatory data into segments and computing the spectrum 
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associated with each segment. The plot essentially shows how the spectral properties 
of the data change with time. 

Using standard Fourier Transform analysis, the spectra associated with each 
time segment have a frequency resolution (Af) which is approximately the inverse of the 
segment duration. 

The model used in this embodiment has two modes. The natural frequency of 
mode 1 changes linearly with time from 90 Hz at the start of the simulation (t = 0) to 1 10 
Hz at the end (t = 100s). The Q level of this mode is 100 which implies a modal 
bandwidth of approximately 1 Hz. 

The digital computer simulation runs over a period of 100 seconds with a time 
step of 0.0002 seconds. The response data seen is typical of what might be measured 
from a vibration survey of an engine in which the engine speed is varying. 

Mode 2 has a natural frequency which is equal to that of mode 1 , plus an 
oscillatory component. The Q level of this mode is 20. The respective modal root- 
mean-square (rms) response levels are approximately the same. Figure 10 illustrates 
these characteristics. 

The objective is to identify the modal bandwidth and damping associated with 
mode 1. 

In this case the zmod data may be analysed to give an estimate of the natural 
frequency profile of mode 1, for example by a standard modal curve fitting approach 
and/or manually cursoring likely estimates prior to conducting a least-squares type fit to 
the data to determine a smooth frequency profile. This profile is shown in Figure 1 1 
and is equal to fn(t) = a)n(t)/2TT. The shift frequency is chosen to be 50 Hz; i.e. Ws = 
2Tr.50 and t* is chosen to be zero. 

The equations described in the general theory section above are then applied to 
yield a carrier signal with variable frequency, in this embodiment using the upper 
sideband formulation. The resulting carrier frequency is shown in Figure 12. The 
carrier phase is shown in Figure 13. Figures 14 and 15 show the carrier signal over an 
interval of 0.04 seconds at the start and end of the simulation respectively. 

The response data Is bandpass filtered between 80 and 120 Hz and 
subsequently used to modulate the variable frequency carrier signal computed above. 
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To illustrate this process, various expanded plots of the processed signal are 
shown over the time interval of 50.0 to 50.3 seconds in Figures 16-19. Figure 16 shows 
the raw response. Figure 17 shows the response after it has been bandpass filtered 
between 80 and 120 Hz. Figure 18 shows the carrier signal. Figure 19 shows the 
modulated signal resulting from multiplying the bandpass filtered response (Figure 17) 
and the carrier signal (Figure 18). 

The zmod plot of the resulting modulated signal is shown in Figure 20, showing 
both upper and lower sideband data. As can be seen in Figure 20, the upper sideband 
data now appears stationary with a constant frequency of 140 Hz in so far as mode 1 is 
concerned. The mode 2 still has a natural frequency which oscillates around that of 
mode 1 . 

If mode 2 had had a natural frequency profile which was simply offset by a 
constant amount from that of mode 1 , then the transformed data would have effectively 
made stationary both modes. 

The transformed data may now be processed using standard techniques as 
discussed in relation to the first embodiment above. For example, the PSD may be 
computed as shown in Figure 21 . Further analysis of this PSD demonstrated a correct 
estimation of the modal damping implicit within the model being simulated. 

If the excitation force is also known (or measured), as is the case in this 
example, then it too may be processed by the method of the present invention to 
generate a modulated signal of the bandpass filtered force. A standard cross-spectral 
based analysis between the respectively modulated force and response signals will 
yield frequency response function data around the mode which has been made 
stationary and whose modal frequency now appears at 140 Hz (90Hz initial natural 
frequency of Mode 1 + 50Hz shift frequency). This frequency response function data 
contains both magnitude and phase data which aids the accurate determination of the 
modal bandwidth and damping. The magnitude and phase of the frequency response 
function resulting from processing this data are shown in Figures 22 and 23 
respectively. 

Third Embodiment 

A third embodiment of the present invention is described in relation to an engine 
order excitation. Again the upper sideband formulation is used here. 
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A digital simulation is made of an engine order (EO) traversing a mode whose 
natural frequency is also varying with time. The zmod or waterfall associated with the 
resulting response data is shown schematically in Figure 24. The variation of the modal 
frequency is also shown to illustrate what would be seen in reality due to low level 
random (asynchronous) excitation superimposed on the engine order excitation. The 
modal frequency varies linearly over the period of interest from 90 to 110 Hz. The 
model mode bandwidth is 1 Hz (i.e. associated Q = 100). The EO frequency varies 
linearly from 80 to 120 Hz over the same time interval. 

The forced response due to the engine order is determined by some form of 
tracking which determines the response at or near EO frequencies. The response 
amplitude profile for this simulation is shown in Figure 25. Analysis of this response 
profile to determine the -3 dB modal bandwidth (related to the damping) yields a value 
of 1 .9 Hz which is almost twice the true value. This error arises due to the natural 
frequency variation with time, which causes the response profile to distort relative to 
what would have resulted had the modal frequency remained constant. 

From the zmod plot, the natural frequency time profile may be estimated and a 
variable frequency carrier signal determined as described above. The response is 
bandpass filtered around the EO frequencies and the resulting signal is modulated with 
the carrier signal as before. The resulting modulated signal is bandpass filtered around 
the upper sideband frequencies of the modulated signal. 

In this embodiment, the shift frequency is chosen to be 50 Hz, with the upper 
sideband frequencies associated with the EO frequencies varying linearly from 140 to 
160 Hz. The tracked response profile of this response is shown in Figure 26 and may 
now be analysed. Analysis of this response profile gives the correct modal bandwidth 
of 1 Hz. The frequency of the tracked response of the upper sideband of the 
transformed data is shown in Figure 27 which demonstrates the correct and continuous 
nature of the transformation. 



